Large spiral and target waves: Turbulent diffusion boosts scales of pattern formation 
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In absence of advection, reaction-diffusion systems are able to organize into spatiotemporal pat- 
terns, in particular spiral and target waves. Whenever advection is present and can be parameterised 
in terms of effective or turbulent diffusion D, , these patterns should be attainable on much greater, 
boosted lengthscale. However, so far, experimental evidence of these boosted patterns in turbulent 
flow was lacking. Here, we report the first experimental observation of boosted target and spiral pat- 
terns in an excitable chemical reaction in a quasi two-dimensional turbulent flow. The wave patterns 
observed are ~ 50 times larger than in the case of molecular diffusion only. We vary the turbulent 
diffusion coefficient D* of the flow and find that the fundamental Fisher-Kolmogorov-Petrovsky- 
Piskunov (FKPP) equation Vf oc \fD* for the asymptotic speed of a reactive wave remains valid. 
However, not all measures of the boosted wave scale with D, as expected from molecular diffusion, 
since the wavefronts turn out to be highly filamentous. 



Pattern formation in reaction-diffusion-advection 
(RDA) systems is an important process in many natu- 
ral and man-made systems, e.g., plankton growth and 
iron fertilization in the ocean pQ, dispersion of pollutants 
in the atmosphere, and optimal mixing in chemical re- 
actors P] . Spiral and target waves have been observed 
on small scales in various active media, e.g. in chicken 
retina [F , cardiac tissue |1 or chemical reactions (5J [5] . 
From a geophysical viewpoint it is of crucial interest 
if these reaction-diffusion patterns can also be found in 
large scale systems involving turbulent advection, as for 
example plankton dynamics in the ocean affecting CO2 
absorption |TJ [51 [7] . Theoretically, the appearance of spi- 
ral and target waves should be possible in RDA systems 
whenever the advection term can be parameterised as a 
global diffusion coefficient |B| . However, so far, experi- 
mental evidence of these patterns in turbulent flows is 
lacking. Despite the importance of pattern formation in 
RDA systems only very few laboratory experiments on 
turbulent fluid flow involve reaction kinetics (9] , and to 
our knowledge, none has considered excitable kinetics so 
far. Considerable numerical and experimental effort has 
focused on cellular and chaotic flows due to the simpler 
realization |2l ITUHl"2~] . In this Letter, we show experimen- 
tally that pattern formation, in particular, spiral and tar- 
get waves can occur in turbulent fluid flows and we find 
that the front expansion is limited by the FKPP equa- 
tion. 

We create a quasi two-dimensional turbulent flow using 
the Faraday experiment |141I15| . i.e. we vertically vibrate 
a circular container of 30 cm diameter filled with 2 mm of 
an excitable cyclohexandione and ferroin based Belousov- 
Zhabotinsky reaction (BZ) [IB] (see methods summary 
and supplementary Fig. SI [IB]). The dynamics of this 
chemical reaction can be well observed in the visible 
range due to the oxidation of the reddish catalyst ferroin 
[F e(phen)\ +1 \ to the blue ferriin [Fe(phen)l + ] [17. We 
vary the intensity of the turbulence and thus the turbu- 
lent diffusion constant |12| by altering the amplitude 



ao of the acceleration and the frequency / of the vertical 
forcing. 

Figure[T] shows examples of the boosted patterns in the 
turbulent flow. The upper panel, 1-3, shows an image 
sequence of a spontaneous boosted spiral and the lower 
panel, 1-3, a spontaneous boosted target wave (supple- 
mentary movies Ml and M2 [13!). Without any fluid flow 
the much smaller usual target and spiral patterns can be 
observed which are shown for comparison on the right 
(image 4). The boosted patterns are a very robust phe- 
nomenon and were found for a large range of forcing pa- 
rameters, / = 30-140 Hz, a,Q = 0.6-2.5 g, g being the 
gravitational constant. The temporal persistence of the 
target patterns varies from some minutes for high forc- 
ing amplitudes ao, to up to one hour for lower ones. The 
probability for a target to form is higher for lower forcing. 
This is most likely related to more long-lived structures 
in the fluid flow [TS] that favor the occurrence of a per- 
turbation that is persistent and big enough to trigger a 
new wave |19| . Usually, but not always, target waves are 
triggered at the border of the container. Spiral waves 
form spontaneously, most often created by the breakup 
of target waves due to interactions with the turbulent 
fluid flow or the boundary, but they can also be created 
intentionally by an abrupt short interruption of the forc- 
ing. Figure[l] (upper panel, image 2) shows the trajectory 
of a spiral tip in time. The temporal persistence of the 
spiral is limited due to the complex movement of the tip 
[2"U] since it eventually hits the border or another pattern, 
causing the spiral to vanish (supplementary movies Ml 
and M3[13]). Qualitative observations suggest that the 
displacement of the spiral tip is a superposition of a ran- 
dom movement due to the filamentary structure of the 
front and a migration along the border of the container 

For a quantitative analysis of the periods of the 
boosted spirals we varied the turbulent diffusion of the 
flow. This was achieved by changing only the forc- 
ing amplitude qq leaving the forcing frequency, and 
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FIG. 1. (Color online) Boosted spiral and target patterns in a turbulent flow. Greyscale indicates concentration of ferriin, 
Fe(phen)^ + . Upper panel, 1-3: Image sequence of boosted spiral with At » 3.6 s, / = 70 Hz, ao ~ 1.8 g, [H2SO4] = 1.2 M, 
period of spiral T = (13 ± 1) s. 2nd image: Trajectory of spiral tip in time (color code: early position blue, late position red). 
Lower panel, 1-3: Image sequence of boosted target wave with At « 12.4 s, / = 50 Hz, ao ~ 1.2 g, [H2SO4] = 0.6 M. Both 
patterns form spontaneously and are persistent phenomena that can last from a few minutes up to one hour. For corresponding 
movies (Ml, M2) see supplementary material |13|. Right: Three close-ups show molecular-diffusion- induced spiral and target 
patterns in absence of fluid flow in the same container. Note the large difference in scales between these usual and the boosted 
patterns. 



thus the Faraday wavelength A f constant [141 |2"2"| (/ = 
50 Hz, [H2SO4] = 0.6 M, supplementary example movie 
M3 [13]). The periods of the boosted spirals at / = 50 Hz 
are in the range T = 30-50 s for all forcing amplitudes 
with a slight tendency towards higher periods for stronger 
forcings. This might be explained by the augmentation 
of the width of the boosted autowaves such that the spi- 
rals seem to be restricted by their own tail [20 . This 
self-restriction could also explain why the period of the 
molecular-diffusion-induced spiral, T mo i = 18-25 s, was 
somewhat lower. Further, in order to prevent the spi- 
ral to drift, we pinned its tip to a round obstacle of 
54 mm diameter, placed in the middle of the container. 
These pinned spirals last for up to ~ 1 h (see supple- 
mentary data, movies M4, M5 and M6 [13]). In addition 
to the spiral and target patterns we also observe double 
spirals with two free curling ends (supplementary movie 
M7[I3]), as well as up to 3 simultaneously existing spi- 
rals. All reactive waves had the typical characteristics 
of autowaves, in particular, they annihilate when they 
meet. 

Figure [2](a) and inset (b) show that the FKPP rela- 
tion for the front velocity v f remains valid for well devel- 
oped boosted target waves in the quasi two-dimensional 
turbulent flow, i.e., Vf = D*/T reac , T reac being the 
reaction timescale [S]. Surprisingly, the boosted data 
points agree with the prediction derived from the FKPP 
equation using only measurements from experiments 
with molecular diffusion: The solid line is the solu- 
tion of the FKPP relation, where the typical reaction 



timescale r reac was estimated from the velocity mea- 
surement of the molecular-diffusion-induced target wave 
to be T reac = (0.8 ± 0.3) s and the molecular diffu- 
sion coefficient was estimated from the literature to be 
Dmoi ~ (1.3-2.0) 10" 3 mm 2 /sdH 121 IM|- Theoretically, 
when the reaction timescale is small in comparison to 
the timescale of the fluid flow, the front velocity Vf is 
bounded by the unidirectional root-mean-square velocity 
of the flow instead of obeying the FKPP relation [25 . In- 
set (c) shows that in our experiments this limit is only 
approached for low forcing. We noted that the variation 
of the front velocity is related to the interval in between 
successive waves which suggests that they might obey a 
dispersion relation analogue to usual target waves [26J . 

In Fig.|2](e) the measured turbulent diffusion coef- 
ficient is plotted ctS £t function of the estimated 
Reynolds number for different forcing amplitudes. The 
turbulent diffusion increases approximately linearly with 
the Reynolds number as expected, and mixing is en- 
hanced. At these Reynolds numbers the flow is turbu- 
lent as can be seen in an exemplary energy spectrum 
(Re « 120) revealing a double cascade and a Kolmogorov 
type scaling (oc fc~ 5/3 ) in inset Fig.|](f) [TSJ [27. The 
turbulent diffusion coefficients were estimated from 
measurements of the absolute dispersion Ait) shown in 
Fig. [2] inset (e), by a fit to the regime of linear growth. 

Despite the validity of the FKPP prediction for the 
front speed, Fig. [3] demonstrates that the boosted target 
waves do not entirely behave like their molecular diffusion 
counterparts. An important difference is the complex fil- 
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FIG. 2. (Color online) Front velocity of reactive waves in dependence of turbulent diffusion, (a) The velocity of the target 
wave fronts Vf scales with \JD, and follows the FKPP prediction Vf = 2^/ D/T reac (solid line). The time constant of the 
reaction r reac = (0.8 ± 0.3) s was derived from the molecular case (circle) but adjusts also well for the turbulent data (crosses). 
Dashed lines indicate the error bounds estimated from the standard deviation of the velocity measurements from the molecular- 
diffusion-induced target wave. Inset (b) shows a close up of the turbulent data pairs, (c) Target front velocity Vf vs. turbulent 
root-mean-square velocity in one direction v' = « r ms/v2, both normalized to the front velocity v mo i of the molecular-diffusion- 
induced target wave, (d) The measured diffusion coefficients are shown as a function of the Reynolds number Re = v rras \Fjv 
indicating the turbulence strength, where v is the kinematic viscosity of the fluid. Inset (e) shows exemplary the absolute 
diffusion for the flows with Re w 43, Re ~ 120 and Re « 194 and the linear fit for estimation of the turbulent diffusion 
coefficient. Inset (f) shows an exemplary energy spectra of the flow for Re « 120. A double cascade and a regime with a 
Kolmogorov type scaling (Ek oc k~ 5 ^ 3 ) can be distinguished, hp is the typical Faraday wavenumber. 



amentous structure of the reaction front which is related 
to the small scale stretching and folding processes in the 
turbulent dynamics (Fig.|3j(a), (b) and Fig, [l) P EH HQ. 
For smaller turbulent diffusion (Fig.[3](a)) the filamen- 
tary structure increases due to two distinct processes: 
First, the increase of the length and persistence of the 
filaments can be explained by coherent flow structures, 
i.e. little eddies and jets, that order the flow on timescales 
longer than the reaction time r reac . An imprint of the 
filaments can be seen in the ferriin concentration profiles 
(Fig.[3](c)). The peaks of high concentration ahead of the 
front show the intermittency of the turbulent diffusion 
process on these spatiotemporal scales. For higher tur- 
bulent forcing the fronts are less intermittent (Fig.[3](d)). 
Second, the sharper and more pronounced appearance of 
the filaments can be explained by the Damkohler num- 
ber, Da = Tfi ow /T reac , the ratio of the typical timescales 
of the flow and the reaction. The flow timescales were 
estimated to be the ratio of the Faraday wavelength and 
the root-mean-square flow velocity, Tfi ow — A/ /v rms . Da 
varied from Da 0.4 for the highest forcing to Da ss 1.8 
for the lowest (supplementary Fig. S2 [13 ). For small Da, 
the fluid flow is fast compared to the reaction timescale 
which causes the front to be smoother in agreement to 
what we find for strong forcing. For large Da, and thus 
lower forcing, the front appears sharper and its velocity 



approaches the root mean square velocity in one direc- 
tion, v' = v rms /y/2 |25], [28] . This limit is reached in our 
experiments for small forcings as is reflected by inset (c) 
in Fig.[| 

Furthermore, in Fig.[3](a),(b), it is easy to observe by 
eye the differences in the target front velocities, the fre- 
quencies of spontaneous target formation and the target 
widths for the two extreme cases of the measured turbu- 
lent diffusion. In order to quantify the dependence of the 
width on the turbulent diffusion, Fig.[3](g) depicts the 
mean profiles of the boosted target for the two turbulent 
diffusion coefficients. These measurements were repeated 
for all turbulent diffusion coefficients (Fig.[3](h)). While 
the full width u>2 of the boosted target waves increases ac- 
cording to W2 oc y/D* , as expected for an ideal reaction- 
diffusion system [29], the width of the rising edge w\ does 
not change within the error of the measurement. A pos- 
sible explanation for this unexpected behavior of Wi is 
the intermittency of the mixing process: Averaging over 
many sharply defined filaments could give a similar width 
for the mean profile as the average over a smoother and 
broader front. This indicates that for low forcings and 
on the timescales of the fast forward reaction occurring 
at the leading edge of the front mixing might not yet be 
well defined by a diffusive process. According to this pic- 
ture, W2 augments diffusively as the backward reaction 
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FIG. 3. (Color online) Front characteristics of boosted target 
waves, (a) , (b) Space-time plots of boosted targets for D» ~ 
5.4 mm 2 /s (a « 1.3 fl ) and £>* « 30.0mm 2 /s (a « 2.230), 
arrows indicate the direction of front propagation (supple- 
mentary movies M2 and M8 |13j). The target waves are nar- 
rower, slower and more filamentous for the smaller diffusion 
coefficient, (c), (d) Ferriin concentration, [Fe{phen)g + ], along 
a line at three different instances of time, At w 6.4 s, for 
D* w 5.4 mm 2 /s and D, « 30.0 mm 2 /s respectively, (e), (f) 
Ferriin concentrations for the same values of D, at three dif- 
ferent points in space (Ax « 53 mm, Aa; « 80 mm). (g) The 
mean profile of the target waves for both diffusion coefficients 
estimated by averaging over all targets measured, (h) Dif- 
ferent widths wi and U12 of the profile in dependence of the 
diffusion coefficient D, . The full width W2 of the target wave 
grows with \fDl as expected while the width of the rising 
edge wi stays constant. 



at the tail of the front is much slower and sees a well de- 
veloped diffusive process. Timescales of the forward and 
the backward reaction can be estimated as the times of 
rise and fall of the ferriin concentration in Fig.[3](e), (f). 

In summary, we conclude that complex spatiotemporal 
patterns, such as target and spiral waves, occur in turbu- 
lent fluid flows as was shown experimentally. Measuring 



turbulent diffusion coefficients and the reaction front ve- 
locities at various Reynolds numbers we find that they 
obey the FKPP relation for reaction-diffusion systems. 
The overall patterns resemble those of their molecular 
counterparts, however, an important difference is the fil- 
amentary appearance of the front which leads to an un- 
expected scaling of the front width. We suggest that this 
phenomena can be understood by the existence or ab- 
sence of coherent structures in the flow that are known 
to exist in many turbulent flows. We expect our results 
to increase the attention on pattern formation in systems 
where excitable dynamics evolve in turbulent flows, such 
as plankton growth in the ocean where a ring-like struc- 
ture, similar to a target, has been reported |3Tj| I31j. 

This work was supported by the Ministerio de Edu- 
cation y Ciencia under Research Grants No. FIS2010 — 
21023. A.v.K. and F.H. receive funding from FPU, No. 
AP-2009-0713 and AP-2009-3550. 



[3 
[4 
[5 

[6 

[7 

[8 

[9 

[10 
[11 

[12 

[13 

[14 

[15 

[16 

[17 
[18 
[19 
[20 



Corresponding author: alexandra.vonkameke@usc.es 
Smetacek, V. et al, Nature, 487, 313-319 (2012). 
Neufeld, Z. & Hernandez-Garcla, E. Chemical and bio- 
logical processes in fluid flows. Imperial College Press, 
London (2010). 

Yu, Y. et al, Nature, 109, 2585-2589 (2012). 
Luther, S. et al, Nature, 475, 235-241 (2011). 
Kapral, R. & Showalter, K. Chemical waves and patterns. 
Kluwer Academic Publishers, Dordrecht (1995). 
Perez-Munuzuri, V., Aliev, R., Vasiev, B., Perez- Villar, 
V. & Krinsky, V.I., Nature, 353, 740-742 (1991). 
Huhn, F. et al, Geophys. Res. Lett. 39, L06602, (2012). 
Grigoriev, R. & Schuster, H.G. Transport and mixing in 
laminar flows. Wiley- VCH, Berlin (2011). 
Ronney, P.D., Haslam, B.D. & Rhys, N.O., Phys. Rev. 
Lett., 74, 3804 (1995). 

Boehmer, JR. & Solomon, T.H., EPL, 83, 58002 (2008). 
Arratia, RE. & Gollub, JR., Phys. Rev. Lett., 96, 024501 
(2006). 

Majda, A.J. & Kramer PR., Phys. Rep., 314, 237-574 
(1999). 

See Supplemental Material at for methods summary, sup- 
plementary figures and movies. 

Haslam, B.D. & Ronney, P.D., Phys. Fluids, 7, 1931- 
1937 (1995). 

von Kameke, A., Huhn, F., Fernandez-Garcia, G., 
Munuzuri, A. P. & Perez-Munuzuri, Phys. Rev. Lett., 
107, 074502 (2011). 

von Kameke, A., Huhn, F., Fernandez-Garcia, G., 
Munuzuri, A. P. & Perez-Munuzuri, Phys. Rev. E, 81, 
066211 (2010). 

Kurin-Csorgei, K., Zhabotinsky, A.M., Orban, M. & Ep- 
stein I.R., J. Phys. Chem., 100, 5393-5397 (1996). 
Perez-Munuzuri, V. & Huhn, F., Nonlin. Processes Geo- 
phys., 17, 177-186 (2010). 

Foerster, P., Muller, S.C. & Hess, B., Proc. Natl. Acad. 
Sci. USA, 86, 6831-6834 (1989). 

Mikhailov, A.S. Foundations of synergetics I: Distributed 
active systems. Springer- Verlag, Heidelberg (1990). 



•5 



[21] Mikhailov, A.S. & Zykov, In: Kapral, R. & Showalter, 
K. Chemical waves and patterns. Kluwer Academic Pub- 
lishers, Dordrecht (1995). 

[22] Tufillaro, N.B., Ramshankar, R. & Gollub, J.R, Phys. 
Rev. Lett., 62, 422-425 (1989). 

[23] Mori, E., Schreiber, I. & Ross, J., J. Phys. Chem., 95, 
9359-9366 (1991). 

[24] Kuhnert, L., Krug, H.-J. & Pohlmann, L., J. Phys. 
Chem., 89, 2022-2026 (1985). 

[25] Brandenburg, A., Haugen, N.E.L. & Babkovskaia, Phys. 
Rev. E, 83, 016304 (2011). 

[26] Dockery, J.D., Keener, J. P. & Tyson, J. J., Physica D, 



30, 177-191 (1988). 
[27] Boffetta, G. & Ecke, R.E., Annu. Rev. Fluid Mech., 44, 

427-451 (2012). 
[28] Koudella, C.R. & Neufeld, Z., Phys. Rev. E, 70, 026307 

(2004). 

[29] Murray, J.D. Mathematical Biology. Springer- Verlag, 

Heidelberg (1989). 
[30] Dubois, D.M. Mem. Soc. Roy. Sci. Liege, 6th series, VII, 

75-82 (1975). 
[31] Wyatt, T., Mar. Biol.22, 137-158 (1973). 



